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Abstract 

We investigate the origin of order in the low-lying spectra of many-body 
systems with random two-body interactions. Our study based both on an- 
alytical as well as on numerical arguments shows that except for the most 
J-stretched states, the ground states in the higher J-sectors are more orderly 
and develop larger energy gaps than the ones in the J = 0-sector. Due to 
different characteristic energy scales in different J-sectors the J = ground 
states may predominate only when all the states are taken together. 
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I. INTRODUTION 



Nature of the mechanism generating order out of randomness constitutes one of the 
most fundamental issues of the contemporary physics. Theories based on various versions of 
ensembles of the random matrices provide one possible theoretical frame for studying such 
effects. The Gaussian orthogonal ensemble (GOE) constitutes the most common reference. 
The related concepts originate from nuclear physics and prove very useful in the area 
of strongly interacting Fermi systems or in quantum chaos 0,^]. At present they finds even 
broader applications in such diverse fields like the brain research 0, econophysics and 
most recently in the "Real- World" networks or graphs Utility of the standard random 
matrix theory (RMT) results form the fact that a potential agreement reflects the generic 
properties of a system and thus in many cases it provides an appropriate null hypothesis. 
From this perspective the deviations are even more interesting as they can be used to quantify 
some properties which are nonrandom and thus system specific. 



In this context the recently identified |[T3JTl[| preponderance of the J = ground states 
in strongly interacting Fermi systems, such as atomic nuclei, arising from random two- 
body interactions seems to indicate the effect reflecting a 'sparser connectivity' than just 
pure random. Several closely related issues have also been addressed in the context of 
mesoscopic []T2[ and randomly interacting many-spin systems |^^. One purpose of the 
present investigation is to identify the origin of the related physically relevant deviations 
from standard RMT and to quantify their character. Since it was nuclear physics which 
gave birth to RMT we believe that the present example, even though addressed in the 
nuclear context, may also stimulate much broader activity and understanding of similar 
effects in other areas. 
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II. STATISTICS OF MATRIX ELEMENTS 



Our theoretical framework is thus analogous to this of ref. |10]. Then schematically, 
indicating nevertheless all the relevant ingredients, the interaction matrix elements v^^, of 



good total angular momentum J in the shell-model basis \a) can be expressed as follows WA 



V. 



J' ii' 

The summation runs over all combinations of the two-particle states \i) coupled to the 
angular momentum J' and connected by the two-body interaction g. g^-, denote the radial 
parts of the corresponding two-body matrix elements while Cj""' globally represent elements 
of the angular momentum recoupling geometry. Structures analogous to eq. (|l|) can be 



identified in various other areas. The quantum open systems |I6[ or the neural nets |T5 
provide immediate examples. 

In statistical ensembles of matrices the crucial factor determining the structure of eigen- 
spectrum is the probability distribution -Py(f ) of matrix elements [|17|. Especially relevant 
are the tails of such distributions since they prescribe the probability of appearance of the 
large matrix elements. From the point of view of the mechanism producing the energy 
gaps they are most effective in generating a local reduction of dimensionality responsible for 
such effects. In principle, the probability distribution of the shell model matrix elements is 
prescribed by their general structure expressed by the eq. (0), provided the probability dis- 
tributions of both g^-i and Cj"^' are known. In general terms this structure can be considered 
to be of the form 



N 



v=j:v, (2) 

2=1 

and each to be a product of another two variables denoted as Cj and Gi. By making use 
of the convolution theorem [0 the probability distribution -Pv'(f) that V assumes a value 
equal to v can be expressed as: 

Pviv) = F-'[FiPvAvi)) ■ F{Pv,{v2)) ■ ... ■ F(Pv.K))], (3) 



where F denotes a Fourier transform, F ^ its inverse and PyXvi) the probabihty distributions 
of individual terms. Taking in addition into account the fact that 



PvAv^)= [ dg.PGMPcA-)T^ (4) 



one can exphcitely derive the form of Pv{v) in several cases. Assuming for instance that 
all the above constituents are identically Gaussian distributed (then, according to eq. @, 



Pvi{vi) = Ko{\vi\)/iT and thus F(Py^(t>i)) = I/vT+cIj^ ) one arrives at 

^v[0) 2(^-i)/2r(Ar/2)0F ' 
where K stands for the modified Bessel function. Asymptotically, for large v, this leads to 

Py(w) ~ exp(— |t>|) |t>|^^^~\ (6) 

For such a global estimate the identical Gaussian distribution of gf-, is consistent both 
with the Two-Body Random Ensemble (TBRE) [|l^] and with the Random Quasiparticle 
Ensemble (RQE) ||10|. The only anticipated difference originates from the fact that in 
the second case the variance of the distribution drops down with J' like the inverse of 
2J' + 1 which is expected to result in a smaller effective as compared to TBRE. By 
contrast, in both versions of the above random ensembles the geometry expressed by Cj""' 
enters explicitely. However, the complicated quasi-random coupling of individual spins is 
believed I^Oj to result in the so-called geometric chaoticity 0. For the extreme values of 
J the underlying selection rules may however impose severe constraints in achieving such a 
limit. Below we therefore explicitly verify its range of applicability. 

III. THE MODEL AND RESULTS 



The model to be quantitatively explored here consists, similarly as in ref. |T0[, of 6 
identical particles (all single-particle energies are set to zero) operating in the sd shell. 
Fig. 1 shows distributions of the corresponding geometrical factors Cj°"' for a a' and for 
several relevant values of J and J'. 



As one can see, the Gaussian may be considered quite a reasonable representation of 
the distribution of such factors for all combinations of J and J' shown, with one exception, 
for those which involve J = 0. In this later case the distribution of Cj?j"' resembles more a 
uniform distribution over a finite interval located symmetrically with respect to zero. One 
principal reason for this fact is that the 6j symbols which enter Cj""' are here more selective. 
These empirical facts justify well the estimates of Pyiv) based on eq. (|]) for J 7^ and not so 
well for J = 0. More appropriate in this particular case is to assume a uniform distribution 
of Cj?j"' over an interval confined by say — cq and cq, i.e., Pc'.(q) = l/2co, retaining PgXQi) 
in its original Gaussian form of course. By making use of eqs. @ and @) one then obtains 



Pv{v) = - 

vr Jo 



00 r 



erf (cqcu/v^) 



CqUJ 



N 



cos{LJv)du (7) 



which for large v behaves like 

Pv{v) ~ exp(-|tf ). (8) 

An explicit calculation of the distribution of the shell model off-diagonal matrix elements 
for the various J-values based on the present model with two-body matrix elements drawn 
from RQE (TBRE results in similar relations among different J-sectors though the distribu- 
tions are somewhat broader as compared to RQE) confirms the above analytical estimates 
as is illustrated in Fig. 2. 

Indeed, such a distribution in the J = sector resembles more a Gaussian and the large 
V tails of this distribution drop down faster as compared to the J 7^ 0-sectors where this 
asymptotics is exponential (eq. (P)). At the same time the J 7^ sectors are dominated by 
very small matrix elements to a larger degree than J = 0. The probability of appearance of 
a large off-diagonal matrix element which in magnitude overwhelms the remaining ones is 
thus greater for J 7^ than for J = 0. Such an effective reduction of the rank in the former 



case is expected to result in a stronger tendency to localization as compared to GOE |[r?| , |21 



The corresponding characteristics can be quantified in terms of the information entropy 

Mj 

^/ = -EKrinKr (9) 



of an eigenstate labelled by / from the J-sector. The coefficients af„ denote the eigenvector 



components in the basis \a >. Such a mean field basis offers an appropriate reference ||22 
for the present purpose. Since the definition of involves the total number of states Mj 
which differ for different J's, before relating the result to the GOE we normalise K/ to the 
GOE limit of this quantity |2B 



K'aoE = + 1) - ^(3/2), (10) 

where ip is the digamma function. Within our model the so-calculated and RQE ensemble 
averaged quantity for all the states versus their corresponding energies E/ is illustrated in 
Fig. 3. As anticipated, it is not J = whose lowest eigenstate comes out most localised, i.e., 
most regular. The lowest states for several higher J values (like 2 and especially 4) deviate 
much more from GOE. This thus indicates more favorable conditions for the emergence of 
energy gaps for larger J than for J = 0. 

Fig. 3 provides one more information which turns out helpful to properly interpret the 
results. The J = states are spread over the broadest energy interval even though the 
number of states (Mq = 14) is here significantly smaller than for several larger J values 
(Ml = 19, M2 = 33, M3 = 29, M4 = 26). As a result, the average level spacing is a factor of 
few larger for J = than for the remaining ones. 

In Fig. 4 (dashed line) we therefore show distributions of the ground state (£"/) gaps 

s' = {Ei-E',)/D', (11) 



where similarly as in ref. |10 



D'=<Ei-Ei>, (12) 

though here for each J individually. 

As it is clearly seen the J = 0-sector does not significantly distinguish from the remaining 
ones. In view of our investigations presented above one would however expect a reduced 
probability for occurance of the large ground state energy gaps in this particular sector. As 
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the solid lines in Fig. 4 indicate such an effect does indeed take place when D'^ in eq. ([TlD 
is replaced by 

=< Ei,^ - Ei > /{Mj - 2). (13) 

In fact, it seems more appropriate and more consistent with the above global considerations 
to relate the ground state energy gap just to the average global level spacing among the 
remaining states, characteristic for a given J, as expressed by eq. ([T3|). 

Finally one may ask a question why this tendency does not extend to the highest J- 
values. In this connection one has to remember that the off-diagonal matrix elements is 
not the only relevant element. These are the diagonal matrix elements which constitute 
the driving term. Irrespective of the value of J their distribution is always Gaussian-like. 
This can be observed numerically and is consistent with arguments formulated in terms of 
eqs. (|I] - 1^) since the geometrical factors Cj""' entering the diagonal matrix elements are 
always nonnegative. As it is shown in Fig. 5 (dashed lines), increasing however J beyond 4 
results in a significant reduction of the variance of Pv{v) for the diagonal matrix elements 
and consequently a larger fraction of the off-diagonal matrix elements becomes effective in 
mixing the basis states. In addition, due to a smaller number of terms entering the eq. 
for the stretched high- J states the effect of geometric chaoticity is no longer effective and in 
this respect the conditions become similar to those for J = 0. As a result, the distribution 
of off-diagonal matrix elements converts back towards more Gaussian-like shaped, i.e. Ne/f 
becomes larger (for J = 6 not shown in Fig. 2 Neff = 2.67). Superposing the above two 
effects one thus obtains even smaller gaps and even more delocalized states at the edges 
of the spectra as compared to J = 0. In fact the spectral density (solid lines in Fig. (5)) 
becomes even somewhat closer to semicircular in this case. 



IV. SUMMARY 



The present investigation based both on theoretical as well as on numerical arguments 
clearly shows that the many-body problems described in terms of various variants of the 
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two-body random ensembles (like RQE or TBRE) develop quantitatively well identified de- 
viations from the GOE. These deviations can be linked to differences in the distribution of 
matrix elements and quantified in terms of the localisation or of energy gaps in eigenspectra. 
Contrary to the common belief they point to the intermediate total angular momenta as 
those J-sectors whose ground states are ordered most. From this perspective, a predomi- 
nance of the J = ground states [0 can be viewed as a result of mixing states with different 
characteristic energy scales from different J-sectors. It seems also appropriate to notice here 
that the arguments formulated in terms of eqs. provide a more adequate approach 

towards understanding the distribution of matrix elements in realistic nuclear shell-model 
calculations than the ones based on multipole expansion [Q]. Finally, similar parallels be- 
tween the distribution of matrix elements and the structure of eigenspectra relative to the 
GOE can be set in the recent econophysics applications. 

We acknowledge useful discussions with J. Okolowicz, M. Ploszajczak and I. Rotter at the 
early stage of this development. This work was partly supported by KBN Grant No. 2 P03B 
097 16 and by the German-Polish DLR scientific exchange program, grant No. POL-028-98. 
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FIGURE CAPTIONS 



Fig. 1 The normalised distribution of geometrical factors Cj""' entering the off-diagonal 
matrix elements (eq. (0)) for the model of 6 particles in the sd-shell. 

Fig. 2 The probability distributions of nonzero many-body off-diagonal matrix elements in 
different J-sectors drawn from one thousand of RQE samples of two-body matrix elements. 
The energy scale is set by where wji = v'^/{2J' + 1) and wj' determines the RQE mean 
square variance. These distributions are fitted (solid lines) in terms of eq. with treated 
as a fitting parameter. The corresponding best A^'s {N^ff) for each J are listed. By increasing 
the distribution prescribed by eq. quickly approaches (as a consequence of the central 
limit theorem) the Gaussian distribution. In this way the J = distribution is demonstrated 
to be much closer to the Gaussian than the remaining ones whose asymptotic behaviour is 
consistent with a slower, exponential fall-off. 

Fig. 3 The information entropy normalised to its GOE limit {Kf /Kqq^) for all the states 
I from various J-sectors (all positive parity) versus energies {E/) of those states. All the 
quantities are ensemble averaged. The energy units are the same as in Fig. 2. 
Fig. 4 Distributions of ground state energy gaps s"^ as defined by the eq. (|ll|) for successive 
J's. The dashed line uses D'^ defined by eq. ( |12|) while the solid line the one defined by 
eq. (IllD. 

Fig. 5 The probability distributions of the diagonal matrix elements in different J-sectors 
drawn from one thousand of RQE samples of two-body matrix elements (dashed lines) and 
the corresponding spectral densities (solid lines). The energy units are the same as in Fig. 2. 
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